HEP-TH/9509147, UPR-661T 



0^ 



0^ 
(N 



NON-TRIVIAL VACUA IN HIGHER-DERIVATIVE GRAVITATION 
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Abstract. A discussion of an extended class of higher-derivative classical theories of grav- 
' ity is presented. A procedure is given for exhibiting the new propagating degrees of free- 

dom, at the full non-linear level, by transforming the higher-derivative action to a canonical 
second-order form. For general fourth-order theories, described by actions which are gen- 
^ , eral functions of the scalar curvature, the Ricci tensor and the full Riemann tensor, it is 

^ ' shown that the higher-derivative theories may have multiple stable vacua. The vacua are 

, shown to be, in general, non-trivial, corresponding to deSitter or anti-deSitter solutions of 

the original theory. It is also shown that around any vacuum the elementary excitations 
remain the massless graviton, a massive scalar field and a massive ghost-like spin-two field. 
The analysis is extended to actions which are arbitrary functions of terms of the form V^'^i?, 
^T) ', and it is shown that such theories also have a non-trivial vacuum structure. 

> 

^ ! 1. Introduction 

O , 

lO ! In a previous paper [1], we presented a method for reducing a general quadratic theory 

of gravitation to a canonical second-order form. The quadratic action provides an example 
of a higher-derivative theory, where the gravitational equations of motion are higher than 
second order. This means that the theory contains more degrees of freedom than just the 
^ ■ simple massless graviton. By rewriting the action, we showed that quadratic gravitation is 
classically equivalent to Einstein's gravity coupled to a massive real scalar field and a massive 
^ ■ symmetric tensor field describing a spin-two field, with a specific Lagrangian. 
d ' We found that, even for simple quadratic gravity, the reduced action gave highly non-trivial 

potential energy and kinetic coupling terms. This suggests that higher-derivative theories 
have a great deal of structure, not immediately apparent from a simple linear analysis where 
one expands around fiat space. One particular feature, which is the main thesis of this 
paper, is that the transformed theory may have a complicated vacuum structure. Here, by a 
vacuum solution, we mean a stable solution of the second-order theory. In this paper, we will 
only consider vacua in which the auxiliary fields are covariantly constant. It is important 
to note that only once the vacuum of the transformed theory has been identified can the 
nature of the elementary field excitations be discussed. Thus, for instance, the exact mass 
and couplings of the excitations, as well as which excitations are potentially ghost-like, may 
be very different around different vacuum states. In general, they will bear little relation to 
the excitations of the linearized analysis. 
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The vacuum structure of the quadratic theories remained comparatively simple. The only 
stable vacuum was flat space with zero vacuum expectation value for both of the auxiliary 
fields. The elementary excitations were a massive scalar field and a massive ghost-like spin- 
two field. We would like to investigate how other higher-derivative theories can introduce 
a more interesting vacuum structure. A natural generalization is to consider actions which 
are not quadratic but general functions of the curvature tensors. Since these tensors involve 
only second derivatives of the metric, the corresponding field equations can be at most 
fourth-order; that is, the same order as the quadratic actions. We can start by asking how 
many new degrees of freedom we expect in such theories. To give a rough count, we turn 
to the Cauchy problem. First recall that even higher-derivative gravity theories remain 
diffeomorphism invariant. Thus we are always able to transform away four components of 
the symmetric metric, leaving six free components. Since the field equations are fourth-order, 
if the Cauchy problem can be solved we expect to be required to give four initial conditions 
for each independent component of the metric, namely the component itself and the first 
three time derivatives. This would imply that we have at most twelve degrees of freedom. 
Rewriting the theory in a second-order form, since the auxiliary fields are set equal to terms 
involving second derivatives of the metric, we expect that giving the initial conditions of the 
auxiliary fields is equivalent to fixing second and third derivatives of the metric, a total of 
twelve conditions. Thus we find that the auxiliary fields should carry six degrees of freedom, 
just as in the general quadratic theory. This leaves a possible six further degrees of freedom 
in the metric. However, since the metric in the second-order theory obeys Einstein gravity, 
the number of degrees of freedom is reduced to the two helicity states of the usual massless 
graviton. By this rough argument, we expect an action in the form of any general function 
of the curvature tensors to describe, at most, the propagation of a massless graviton plus 
an additional six degrees of freedom; that is, the same as for the special case of the general 
quadratic theory discussed in our previous paper [1]. Be this as it may, the structure of such 
theories is potentially much richer than in the quadratic case. As we will see, the vacuum 
structure is now non-trivial. 

A second generalization is to consider gravity actions which have higher derivatives act- 
ing on the curvature tensors. Since these theories involve more than two derivatives on 
the metric, the corresponding field equations will be generically higher than fourth-order. 
Consequently, we now expect new degrees of freedom in addition to those of the general 
fourth-order theory. Again, we will show that the vacuum structure of such theories is, in 
general, non-trivial. 

Some discussion of reducing both types of generalized theories to a second-order form 
already exists in the literature. That the equations of motion following from actions given by 
a general function f{R) are equivalent to the equations of motion of a scalar field dilatonically 
coupled to gravity, was first shown by Teyssandier and Tourrenc [2] . The equivalence at the 
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level of the action was given by Magnano et al. [3-5], who also rewrote the reduce theory in 
canonical form. This latter group and Jakubiec and Kijowski [6] also rewrote actions given 
by a general function f{R^u) in second-order form, though without the canonical separation 
of the new degrees of freedom. Actions including derivatives of the scalar curvature were 
considered in the context of inflation by Gottlober et al. [7], who again showed, in some 
special cases, the equivalence of the field equations to those describing scalar fields coupled 
to gravity. Still at the level of the field equations, this work was extended to a general class 
of actions by Schmidt [8] and Wands [9]. Possible vacua of the f{R) theory are discussed, 
in terms of the original higher-derivative field equations, by Barrow and Ottewill [10], and 
later in terms of the second-order field equations by Barrow and Cotsakis [11]. 

This paper is organized as follows. In the next three section we will consider actions 
which are general functions first of the curvature scalar only, then of the Ricci tensor, and 
finally of the full Riemann tensor. Rewriting these theories in a canonical second-order form, 
we will find a rich vacuum structure. This will require a careful discussion of the different 
"branches" of the theory, a concept we define below. We will show that the vacuum states are 
in a one-to-one correspondence with the stable, constant curvature, deSitter or anti-deSitter 
solutions of the higher-derivative theory. Furthermore, we will show that if the vacuum of 
the transformed theory has zero cosmological constant, then the corresponding spacetime 
for the original higher-derivative theory is also flat. We will then discuss the elementary 
excitations around a given vacuum and argue that they remain a graviton, a scalar field and 
a ghost-like spin-two field, with masses and couplings depending on the particular vacuum 
and the functional form of the original action. In Section 5 we will consider a class of gravity 
actions that are general functions of the curvature scalar R and the derivatives V'^^R, with k 
a positive integer. We show how to rewrite such theories in second-order form and show that 
they have, in general, a non-trivial vacuum structure. In particular, we present an example 
with two new scalar degrees of freedom, neither of which is ghost-like, coupled to Einstein 
gravity with a stable anti-deSitter space as its vacuum. We briefly present our conclusions 
in Section 6. 

Throughout the paper our conventions are to use a metric of signature ( — h ++) and 
deflne the Ricci tensor as R^, = dxT\, - d^V\^ + V\Tp^^ - T\^T\^. 

2. Actions Given by General Functions of the Scalar Curvature 

As in the case of quadratic gravity, the higher- derivative theory described by actions 
which are general functions of the curvature scalar is classically completely equivalent to the 
canonical second-order theory of a scalar fleld coupled to gravity. This equivalence was flrst 
shown at the level of the action by Magnano et al. [3-5]. We shall derive this equivalent 
theory in a slightly different form, by flrst introducing an auxiliary fleld to reduce the action 
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to second order and then making a suitable conformal transformation. We will point out the 
need to consider different branches of the theory when making the reduction. 
We introduce the auxiliary field in two steps. First we write 

^ 2^ / ^'"^^ [^'^^^ {R-X) + f{X)], (2.1) 

where f'{X) = df /dX. The auxiliary field X has the equation of motion 

r{X){R-X)^0. (2.2) 

Provided we have f"{X) ^ 0, this gives X = R. Substituting back into the action we return 
to the original higher-derivative form. Thus, the reduced action is equivalent to the original 
theory, but only away from the critical points defined by f"{X) = 0. A continuous region of 
X between critical points, where f"{X) never vanishes, will be called a branch. Typically, 
there will be several branches in the theory. Since X gets set equal to R, in terms of the 
original theory, the condition for a critical point can be written f"{R) = 0. Thus, branches 
in X correspond to branches in R in the original higher-derivative theory. 
As a concrete example, let us assume that 

f{R) = R + e-^R\ (2.3) 

where e > 0. Solving f"{X) — imphes that there is a single critical point at X — 0. 
Therefore, the second-order formahsm in terms of the auxiliary field X is valid in each of the 
two branches — oo < X < and < X < oo, but breaks down at the critical point X — 0. 
Since X gets set equal to R, these regions correspond to branches in the space of curvature 
given by cx) < i? < and < i? < oo with a critical point ed. R = 0. 

The second step, which will then allow us to rewrite the action in canonical form, is to 
change variables to a new field A = f'{X). In terms of this new variable, the action becomes 

^ 2^ / ^'"^^ " ^^^^^ + fi^W)]- (2.4) 

Clearly this action is well defined only for those regions where the X action is valid; that 
is, away from the critical points. Therefore, the A formulation is only defined over ranges 
of X for which f"{X) ^ 0. Note that to define the action (2.4) we must be able to invert 
the expression A = f'{X) to give X = X{X). Locally this requires the same non-degenerate 
condition f"{X) ^ 0. Globally, there may still be many different roots when we solve for X 
in terms of A. However, in any given branch of X there is only a single root. Thus there is 
a valid formulation of the theory in terms of A for each branch of X. It is important to note 
that, in each branch, the inverted function X = X{X) which enters (2.4) is different. 
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As a concrete example, consider once again f{R) defined in (2.3). It follows that 

X^l + 3e-^X\ (2.5) 
Note that A > 1. This expression can be inverted to give 

X = ±^v/A^. (2.6) 

Clearly one must take the + root in the < X < oo branch and the — root in the —oo< 
X < branch. This situation is quite generic, as we will see below. 

For a particular reduced theory, corresponding to a given branch of the original higher- 
derivative theory, we can define x — log ^ perform the conformal transformation = 
e^Qui,. Action (2.4) then becomes, 

^ J d'x^[R-l{WxY-e-'^e>^X{e^)-f{X{e^)))\. (2.7) 



2^2 

As promised, we find a theory of Einstein gravity coupled to a scalar field with a particular 
potential dependent on the choice of the original function f{R). As in the quadratic case 
there is one subtlety in defining x- To keep x real, we must have A > 0. When A < we 
must introduce, instead, the field x — log(~'^)) the effect of which is to change the sign 
of the overall normalization of the action as compared with the form given in (2.7) above. 
Further, at the special point A = the action cannot be put in canonical form. In this paper 
we will restrict our attention to the A > case only. 

As a concrete example, consider f{R) defined in (2.3). As pointed out above, in this case 
A > 1 in either branch of X. Therefore, x = ^ogX is always well defined. We point out, 
however, that the range of x is restricted to < x < cxo. 

Note that the scalar field kinetic energy has the usual sign and, hence, x is not a ghost. 
We can now look for the vacua in any given branch of the theory. The vacua of the theory 
described by action (2.7) are defined to be stable solutions of the g^^^ and x equations of 
motion given by 



V„xV„x - \9iM, (Vx) J - ivV(x) (2.8) 

V\ = ^ (2.9) 
dx 

respectively, where the potential is given by 

V{x) = |e-2^ {e^X{e^) - f{X{e^))) , (2.10) 

where we recall that the form of X[e^) depends on the branch in question. As stated 
previously, we will only consider vacua satisfying the covariant constant condition 

V^X = = 0. (2.11) 
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It follows that X is a constant and, from (2.9), that it must cxtremize the potential. Fur- 
thermore, it follows from (2.8) that the vacuum is a space of constant curvature with 

R^6V{x). (2.12) 

Since Q/^i, — e~^g^^ and x is constant, these vacua also correspond to spaces of constant 
curvature with respect to the original metric, with R — e^R. Since the vacuum field x must 
be an extremum of the potential, we have the condition 

e-^^ (2/(X(e^)) - e^X(e>^)) = 0. (2.13) 

We started with a completely general function /, so this condition may generically have 
multiple solutions, including solutions away from x = 0. Further, given a particular solution 
for X, the value of the potential at this point need not be zero; that is, the vacuum generally 
has non-zero cosmological constant. In this case, the curvature scalars R and it! will be 
non- vanishing. 




As a concrete example of all this, we return to the example specified in (2.3). As discussed 
previously, there are two branches where the second order theory is defined; branch A where 
< X < cxo and branch B where — oo < X < 0. The field A, which is given as a function 
of X in (2.5), satisfies A > 1 everywhere and, hence, x = log A is well defined in both 
branches, although its range in restricted to < x < oo. Expression (2.5) was inverted 
to give X as a function of A in (2.6). This expression is branch dependent, being given by 
X = ±(e/ v^) V A — 1 with the positive root chosen in branch A and the negative root chosen 
in branch B. It follows that the form of the potential energy defined in (2.10) also depends 
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Figure 2. Vb{x) for R + e~'^R^ gravity 



on the branch chosen. It is given by 



Vb{x) 



9^/3 



ee 



-2x(eX _ 1)3/2 



(2.14) 



for branch A and B respectively. The two potentials are plotted in Figures 1 and 2. We see 
that each has a stationary point. However, since vacua must be minima of the potential, 
it follows that only branch B has a stable vacuum. This is located at x — which 
is non-zero. Note that VB(ln4) = — e/24, so that the vacuum state has a non-vanishing 
negative cosmological constant. It follows that R — — e/4 and it! = — e, which implies that 
the vacuum state is an anti-de Sitter space both in terms of the metric ^^^^ and the metric 
Qi^,^. Thus we have given an example of a higher-derivative theory which has a new vacuum 
state away from the flat-space solution of ordinary Einstein gravity. It is important to note 
that none of this structure would have been evident if we had made a simple expansion of 
the original action around flat space in terms of hf^u — g^i, — rj/j^i,. In fact, keeping only the 
first non-trivial, quadratic terms in h^,,, we would not even have been aware that there was 
an additional scalar degree of freedom in the theory. The lowest order term in R^ is cubic 
in hfj_u, and so in a quadratic expansion we would only see the usual Einstein term. 

The conclusion is that generic higher- derivative corrections to pure Einstein gravity intro- 
duce completely new vacuum states into the theory. Flat space is no longer a unique point 
in field space, so that if, for instance, we wish to investigate the fundamental excitations in 
the theory, we must start by specifying which vacuum state we are considering. 
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At this point, we would like to present another specific example that further illustrates 
the preceding discussion. Consider the cubic function 

f(R) + Sae-^R^ + e'^R^ (2.15) 

with e > \/3a > 0. We now have 

fix) = 6e-\X + a) (2.16) 

so that the theory has a single critical point at X = —a. It follows that the second-order 
theory is defined on two branches; branch A where — a < X < oo and branch B where 
—oo < X < —a. On either branch we can define a new field A = f'{X), which in this case 
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Note that, since e > \/3a, A is always a positive real number in the range GO 
for both branches. Expression (2.17) can be inverted to give 




Clearly the X+ solution is correct on branch A whereas the X_ solution is to be used in 
branch B. Since A is always positive, the definition of x = In A is valid in both branches. Note 
that X is then restricted to lie in the range ln(l — 3«^e~^) < x < oo. Since the expression 
for X in terms of A, and hence given in (2.18) is different in each branch, it follows from 
(2.10) that the potential energy function also is different in each branch. We have 



and the Va expression is valid in branch A and the Vb expression valid in branch B. The 
two potentials are plotted in Figures 3 and 4. Unlike the previous example, in this case each 
branch contains a single, stable minimum. In branch A, this minimum is located at % = and 
has vanishing cosmological constant V^(0) = 0. It follows that, for this vacuum, R = R = 
and spacetimc is flat with respect to both the 'g^,^ and Qfj^i^ metrics. In branch B, the minimum 
is located at % = ln(4 — 6ae~^) with a non-vanishing negative cosmological constant given 
by V{x) = — e^/6(4e — 6a). For this vacuum R — — e^/(4e — 6a) and, hence, R = — e. It 
follows that the spacetime is an anti de-Sitter space with respect to both metrics g^^, and g^i,. 
Thus we have found an example which, aside from a conventional minimum at R — with 
zero cosmological constant, has an additional minimum with negative cosmological constant. 
Furthermore, there is an unstable maximum with positive cosmological constant. The mass 
of the scalar field is different at each of the different minima. Here a quadratic expansion of 
the action around fiat-space would have identified a new scalar degree of freedom, but would 
never have revealed the presence of a second stable vacuum state. 

It will not have escaped the readers notice that finding the vacua in the second-order 
formalism requires a careful discussion of the branches in the X variable. The branching 
structure was not too difficult in the preceding examples, but as we will show below, it can 
be, and generally is, extremely complicated. However, the following remarks will allow us to 
define a simpler procedure for determining the vacua. Recall that, for a covariantly constant 
scalar x satisfying dV/dx = 0, the corresponding spacetime structure in terms of the metric 
g^iy is that of a space of constant curvature with R = 6e^V{x)- Provided we are not at a 



Va{x) 
Vb{x) 



= ie-V2M^+(x))'(3« + 2^+(x)), 

= |6-V^M^-(x))'(3« + 2X_(x)), 



(2.19) 



where 




(2.20) 
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degenerate point, spaces of constant curvature R and vacuum solutions are in one-to-one 
correspondence, since the relation R = 6e^V{x) not only means that constant x implies 
constant R, but also the reverse. Note that, as a corollary, any vacuum solution with zero 
cosmological constant, namely V{x) = 0, must correspond to flat space in the original higher- 
derivative theory. An important consequence of this result is that, since fiat space is a single 
point in the field space of the original theory, there can be only one vacuum state with zero 
cosmological constant in the second-order theory. We conclude that all vacuum solutions of 
the above type can be found as constant curvature solutions of the equations of motion for 
the higher-derivative g^i, theory. All such vacua can be found directly, without reference to 
any branch structure. Having found a vacuum of interest, one can then proceed in reverse, 
introducing the second-order theory in the appropriate branch containing that vacuum. This 
is often an easier procedure. The gi^i^ equations derived from the higher-derivative action 
(2.1) are 

Rf^ufiR) - y^.uf{R) - (V^V. - g^,V') f\R) = 0. (2.21) 

If we look for solutions of constant curvature, then R^^^ = ^g^^R where i? is a constant. It 
follows that the derivative terms in the equations of motion drop out and we are left with 
the simple condition, first derived by Barrow and Ottewill [10], 

Rf'{R) - 2f{R) = 0. (2.22) 

We see that this is exactly the condition we obtained for stationary points of the potential 
(2.13) aside from a factor of e^"^^. In the latter case, however, the expression was taken 
to be a function of the scalar field x expressed in terms of X, and as such only valid in a 
branch-by-branch sense, whereas here the expression is valid for all curvature R. 

This now provides us with a procedure for finding all the vacua of the theory as well 
as the nature of the excitations around a given vacuum. We start by looking for constant 
curvature solutions of the original higher-derivative equations of motion, solving the equation 
(2.22), which is valid globally. Having identified these points we can then, locally around 
each solution, make a transformation to the reduced theory, choosing, if necessary, the 
appropriate branch oi X{e^). In the transformed frame, the new scalar degree of freedom is 
made exphcit. We can then address questions of local stabihty and the mass of the scalar 
field. The only breakdown of this procedure occurs when the minimum is at a degenerate 
point of the transformation; that is, a point where f"{R) = 0. In this case, any discussion of 
the particle spectrum and local stability must be in terms of the original theory. In general, 
we shall not consider such degenerate points. 

As an example of this procedure, consider the function 



f{R) ^l[R + esm{R/e)] 



(2.23) 
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Where e is a constant. The condition for a degenerate point is then 

f"(R) = -if-i sin (R/e) = 0. (2.24) 

This equation has solutions R — nne where n is any integer. The theory thus has an infinite 
number of branches given by mre < R < {n + l)7re. Rather than treating each branch 
separately, we follow our procedure for finding the vacuum states of the theory by solving 
the condition (2.22) for the constant curvature solutions of the original theory. The condition 
reads 

sin {R/2e) [{R/2e) sin (i?/2e) - cos (i?/2e)] = 0. (2.25) 
The left-hand side of this expression is plotted in figure 5. We have the solutions 

R = 2n7re, 25„e (2.26) 
where n is an integer and 5„ are the solutions of 

5sin5-cos5 = 0. (2.27) 
Explicitly, the first few solutions of (2.27) near 5 = are 

. . . , 5_2^ -6.12, -2.80, 5i K 2.80, 82 ^ 6.12, . . . (2.28) 

Furthermore (|n| — |)7r < \5n\ < \n\TT and as \6\ becomes large the solutions approach 

Sr, = UTT. 



J\ A 




X 27ce 


2536^ 4jc6 


47te^ 25.26 


-27ceV 


/ 25.16 





Figure 5. Condition for a constant curvature solution of | + esin(i?/e)] gravity 

The first observation to make is that the constant curvature solutions at i? = 2n7ie, 
including the point R = 0, all correspond to degenerate points. As such, a second-order form 
of the theory does not exist at these points. Thus, although these arc constant curvature 
solutions, we cannot identify them as vacuum solutions of a canonical second-order theory. 
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However the solutions at i? = 25„e are not at degenerate points. From our previous discussion 
this imphes that they must correspond to vacuum solutions of the canonical second-order 
theory. Nonetheless, we find that each vacuum solution lies in a different branch of the 
theory. 




Figure 6. V{x) for the branch ne < R < 27re 



Let us consider two particular solutions, R = 26ie ^ 5.60e and R — 25_ie = —2Sie ~ 
— 5.60e, in the second-order formulation. For the first solution we must take the branch 
Tie < R < 27ie. Following our usual procedure for writing the theory in a second-order form, 
we first introduce the auxiliary field X and then define 

A = /'(X) = i[l + cos(X/e)]. (2.29) 

However, as usual, we are required to invert this expression to give X as a function of A. It 
follows that 

X = ecos"^ (2A - 1) with Tre < X < 27re. (2.30) 

We note that A is only defined in the range < A < 1. To put the action in canonical 
second-order form we define x = In A. The potential for x (2-10), is then given by 

V{X) = le-'^[X{x)cos{X{x)/e)-esm{X{x)/e)] (2.31) 

where, 

X(x) = ecos~^(2e^-l) with 7re < X < 27re. (2.32) 

We now have that x is restricted to the range — oo < % < 0. The potential is plotted in figure 
6. We find that there is a single maximum of the potential at x ~ —0.12. The value of the 
potential at this point is V'(x) 1.05e, so that R = QV{x) ~ 6.31e and R = 6eW{x) ~ 5.60e. 
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Thus we see that the maximum corresponds to the flat space solution R = 2Sie as expected, 
and we can conclude that this vacuum is unstable. 




Figure 7. V{x) for the branch -27re < R < -ne 



For the solution at R — 2(5_ie, we must take the branch —2ne < R < —ne. The expression 
for the potential is the same as the expression (2.31) above, except that we must now define 

X(x) =ecos-^(2e^-l) with -27re < X < -7re. (2.33) 

The resulting potential is plotted in figure 7. Again x is restricted to lie in the range 
— oo < X < 0. We now find a stable minimum at % ~ —0.12. This gives V{x) ~ — 1.05e, 
R — 6V{x) ~ — 6.31e and R — 6e^V{x) ~ — 5.60e, showing that the minimum of the 
potential does indeed correspond to the constant curvature solution R — 2(5_ie. We can 
concluded that this is a stable vacuum state. In fact, further analysis shows that all the 
vacua it! = 2(5„e with 5„ > are unstable maxima, while those with Sn < are stable 
minima. This final example demonstrates that the vacuum structure of higher-derivative 
theories can be complicated. We find an infinite number of stable and unstable vacua, each 
in a different branch of the theory. 

In summary, we have shown that theories which are general functions of the scalar curva- 
ture R, can be reduced to a canonical second-order form describing Einstein gravity coupled 
to a scalar field with a potential which may have a vacuum structure of arbitrary complexity. 
Further these vacua all correspond to constant curvature solutions of the original theory. The 
advantage of discussing the theory in the second-order formalism is that it is in a canonical 
form where we can easily identify the excitations around the vacuum. The disadvantage 
is that the problem of inverting f'{X) means that we are forced to break the theory into 
regions separated by degenerate points. For this reason it is often easier to identify vacua as 
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constant curvature solutions of the higher-derivative theory and then make a local transfor- 
mation in the region of each solution to investigate the excitations around that vacuum. All 
of the non-trivial vacuum structure arises from non-linear terms in the field equations and so 
is missed in a quadratic expansion of the higher-derivative theory in terms of the fiuctuation 
of the metric around fiat space. Further, any analysis of the particle content of the theory, 
in this case the presence of a new scalar degree of freedom and identification of its mass, can 
only be made once the vacuum is identified. 

3. Actions Given by General Functions of the Ricci Tensor 

In this and the following section we shall extend our discussion to include actions which 
are general functions of first the Ricci tensor and then the Riemann tensor. We shall use the 
same techniques as in the previous section, namely introducing auxiliary fields, showing that 
vacuum states correspond to constant curvature solutions of the higher derivative theory, and 
then investigating the excitations around a given vacuum. The reduction to a second-order 
form was first discussed by Magnano et al. [3-5] and Jakubiec and Kijowski [6] . Here, we shall 
use a slightly different analysis, again pointing out the need for identifying different branches 
of the theory and also giving a canonical separation of the new variables about any given 
vacuum. We shall only give the general formulation without specific examples. Again we 
will find a rich structure of vacua which would be completely missed in a linearized analysis. 
This will re-emphasize the need to identify the vacuum in question before investigating the 
masses and couplings of the elementary excitations of the theory. 

Starting with an action f{R,j,v) we can introduce an auxiliary field and put the theory in 
a second-order form 



-5 = ^1 d'xyTgfiR,,) 



2k 



dXfj^jy 



= d'^V^ K'^ (i?^. - X^,{7T,„)) + /(X^,(7r,,))]. (3.1) 

Again we introduce the auxiliary field in two steps, first writing the action in terms of the 
field X^i, which gets set to it!^,^ on solving its equation of motion. We then define 

n^^" = (3.2) 

and invert the expression to give X^^i, as a function of ir^,^. Note again that the introduction 
of both Xf^i, and tT/^i, requires the non-degeneracy condition 

(jg^ ^ fi-^Kx) _^ Q (3 3) 

dXfj_^dXp^ 
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if the auxiliary field is to be properly eliminated to return to the original action. Furthermore, 
in inverting (3.2), it may be necessary to divide the theory into branches corresponding to 
different possible roots for X^,^. Thus for the variable -k^j^ it may be necessary to introduce 
a collection of auxiliary variable theories each valid for a different branch of R^v. 

We can ask how many degrees of freedom are there in the auxiliary field? We have 
argued that we expect six. If we consider the X^^, formulation for now, we see first that 
we can again derive a spin-two divergence condition. Since the auxiliary field equation of 
motion gives X^^ — R/^i,, we clearly can show, using the Bianchi identity V^G^^ = 0, that 
{X/^i, — jQf^uX^ — 0. Constraining the components of X/j^i,, these four conditions imply 
that we do indeed have six new propagating degrees of freedom. We would like to be able 
to separate these into scalar and spin-two degrees of freedom as was done in the quadratic 
case in [1] . For a general function /, this is not possible since we are not able to obtain the 
necessary trace condition from the g^i, equation of motion. However, as we will show below, 
having identified a vacuum state, we can always make the separation locally in an expansion 
around the vacuum solution. 

What then are the vacua of the f{Rfj,iy) theory? Since we are unable to separate the 
spin-two and scalar degrees of freedom in the auxiliary field, we cannot put the theory in a 
canonical form and look to minimize the potential as we did for the f{R) theory. However, 
we recall that the vacuum states in question arc none other than states where the auxiliary 
field is covariantly constant. This means, since the X^^, equation of motion gives Xj^^, = R^^ 
that they are states of constant tensor curvature. We can also impose the condition that 
the states should have no preferred direction, so that X^^, is proportional to the metric 
Qf^,^. Thus vacuum states are states of constant scalar curvature R^i, — ^Rg^^u with R 
constant. Consequently, one approach to finding vacuum solutions is simply to look for 
constant curvature solutions of the original theory. 

The equations of motion of the original higher derivative theory are given by 

VV;. + V - V,V^/r - V.V.f/ + R,,C + R,.f; - g,.f = 0. (3.4) 

To incorporate the presence of the metric in f{Rfj,u) in deriving the equations of motion, 
we consider / as a function of the mixed index object R^", with all contractions made 
between raised and lowered indices so that the metric does not enter explicitly. Then in the 
above expression f'^" — df/dR^,^, so that, for instance, /^^ = f'J^Qxv Looking for constant 
curvature solutions of the form R^^, — ^g^^R, we obtain the condition 

Rg'{R) - 2g{R) = 0, (3.5) 

where g{R) = f{\R5^^) and g'{R) = dg{R) / dR, and we have used the fact that f'^'^ evaluated 
at the constant curvature solution must be proportional to (5^*^. Note that this condition has 
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exactly the same form as the condition (2.22) we obtained for actions which were general 
functions of the scalar curvature. 

Given a particular vacuum solution satisfying the condition (3.5), we would like to inves- 
tigate the excitations around the vacuum state. The natural way to do this is to consider 
an expansion of the action around a given vacuum R = Rq. Expanding fi^R^y) and keeping 
terms to quadratic order in the curvature only, we have 

I "^^""^ [""o + ai{R- Ro) + |a2,i {R - Rof 



— \Ro + i? + 77— ^-R^ — — -5- {Rfj^iiR^^ — 



6mn^ 



(3.6) 



Here we have evaluated / and its derivative at the vacuum R^^ — \g^i,RQ and defined 



/I 



dR. 



Ro 



(3.7) 



Ro 



using the general symmetric decomposition of the second two expressions. The variables (3, 
mo and m2 are then given by 

p = ai- (02,1 + ^02,2) -Ro, 

_l 

(3.8) 



mo 



1712 



3a2,i + 02,2 ' 
02,2' 



where in writing the final line of (3.6) we have used the constant curvature condition aiRo = 
2ao. We see that we have succeeded in putting the action in a quadratic form, though with 
the addition of a cosmological constant. 

The final line of (3.6) is identical to the general quadratic form we discussed in our previous 
paper [1], except for the addition of a cosmological constant term and a renormalization of 
the gravitational coupling constant by a factor (3. We can thus introduce auxiliary fields x 
and TT^i^ exactly as we did in [1] to give 



2/^2 



d X 



R 



l(Wx) - Imo' {1 - e-^y - '^Roe-'^ 

-G^.Tt'^'' + >2' {n^^un'"' - 7f2)l . (3.9) 
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The only effect of tlie cosmological term is to modify tlie potential for tlie scalar field, adding 
a term ^RqC^^^. We find that the extremum of the potential is now at 3mo^ (e^ — 1) = Rq, 
which, relating the auxiliary field x back to the original curvature, gives R = Rq, as required 
for a consistent expansion. Following exactly the analysis of our previous paper, the auxiliary 
field TTfj^i, satisfies generalized divergence and trace conditions, so it does indeed describe the 
degrees of freedom of a spin-two field. By making a final field redefinition, we can write 
the action in canonical form with explicit kinetic energy terms for the spin-two field. Again 
following the discussion in our previous paper, the spin-two field will have the correct Pauli- 
Fierz limit, but will unfortunately be ghost-like. 

Thus we have shown that locally, around any vacuum of the higher-derivative theory (that 
is a solution of constant curvature), we can expand the theory to identify a new scalar and 
a new spin-two degree of freedom (provided we are not at a degenerate point). The spin- 
two field satisfies divergence and trace conditions as before but importantly, we see that it 
remains ghost-like. Thus generalizing to f{Rni^) actions fails to remove the problem of the 
ghost spin-two degree of freedom. The masses of the degrees of freedom are fixed by the 
form of the function / around the constant curvature solution. 

It is worth noting that we could equally well have done this analysis in the second-order 
form, looking for solutions with constant X^j, proportional to g^iy. Then expanding in small 
X^t, about such solutions, gives a linear coupling between X^j, and R^^^ and quadratic 'mass' 
terms for X^j^. From this form we could then extract the scalar and spin- two parts of X^jy. 
In this sense, the auxihary field X^^, always carries six degrees of freedom, which about any 
given vacuum can be decomposed into a scalar field and a ghost spin-two field, with the form 
of the decomposition changing as we go from vacuum to vacuum. 



4. Actions Given by General Functions of the Riemann Tensor 



Our last generalization is to consider actions which are general functions of the Riemann 
tensor Rx^^^p- As mentioned earlier, the suggestion is that such theories have an additional 
six degrees of freedom, which we know in the quadratic case can be decomposed into scalar 
and spin-two fields. 

We start the discussion by demonstrating that, as in all previous cases, we can introduce 
an auxiliary field to write the action in a second order form. 



1 

2^ 

1 
2^ 

1 
2^ 



d'^x 
d^x 



gf{R\nvp) 

df{Xap^s) 



dX 



X/iup 



{Rxnup — X\ixup) + f{X\nvp) 



J d^xV^ [p^'^^P {Rx^^^p - Xxp,^.p{pap-y5)) + f {Xxp,up{Pap-y5))\ ■ (4.1) 
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Again we introduce the auxiliary field in two stages; first introducing Xxfj,up, which gets set 
equal to the Riemann tensor on solving its equation of motion, and then defining 

In both cases we require the non-degeneracy condition 

in order to be able to eliminate the auxiliary field and return to the original action. When 
using the variable pA^ti/p we may be required to break the theory into branches, introducing 
a collection of second-order theories, each taking a different branch when inverting to find 
Xx^,ua in terms of px^,vp- 

Turning to the vacuum states, we again find that states with constant Xxfj.ucr have constant 
Riemann curvature, since Xx^ua is set equal to Rxi^ua by its equation of motion. Furthermore, 
imposing the condition that there is no preferred direction in spacetime we require that Xx^ua 
is proportional to gi^i^. Therefore, given the symmetries of Xx^ua, we have 

R\lJ.ua — R%ua — ^-^0 {gXuQua — QXaQnu) , (4.4) 

and, hence, we are considering solutions of constant Ricci scalar curvature. 

As before the easiest way to obtain such solutions is from the equations of motion of 
the original higher-derivative action. Again, to circumvent the problem, when deriving the 
equation of motion, of the metric explicitly entering the function /, we consider / as a 
function of R^^" with all contractions made between raised and lowered indices. We then 
derive the equations of motion 

vvv^p.. + vv^/^,^, + \Rr-fi,„r + \Rr^r,p.r - b^uf = o, (4.5) 

where f'/fj^ — df /dRJ^^^ . Restricting to constant curvature solutions of the form (4.4), we 
get the familiar condition 



Rog\R^) - 2g{Ro) = 0, (4.6) 

where now g{Ro) = f^R^xp-vp) ^'^'^ o'i^o) = dgi^Ro) j dRQ, and we have used the fact that, by 
symmetry, /'//(i?A^^p) is proportional to S^.'^S/ -gpyg^". 

To investigate the excitations around a given vacuum we expanding the action about the 
constant curvature solution. We have, keeping terms up to quadratic order in the curvature 
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only, 

+|a2,i {R — Ro)^ + |a2,2 {R^u — \giiuRo) {R'^" — \g^''Ro) 
+1^2,3 {Rx^^up - ^Ro [gxug^^p - gxpg^^A) {R^'""' - ^Ro [g^'g'"' - g^'g'"'])] 

+7 {Rx^.pR'^^'"' - ^R^uR"" + R')] . (4.7) 

Here we evaluate / and its derivatives at the constant curvature solution, defining, for the 
contraction of the derivatives with a tensor Ax/^^p which has the symmetries of the Riemann 
tensor but in otherwise arbitrary, 

/Iro = "0' 

df 



JL 

2^2 



1 2 1 

"o-Ro + R + - — T^R 5- (RnuR^'^ — nR ) 



dRxixiip 



Axp.„p = oiA, 

-Ro 



(4.8) 



A^^A^A.^ar = 02,1 A^ + 02,2 A^.A'^" + a2,3AA^,pA^'^"'', 

-Ro 



(4.9) 



dRr}K\p.dR 1 

with = g^^A^xvp and A = g^'^A^^. Wc have also introduced the parameters 

P 2 2/3 

3a2,i + a2,2 + 02,3 ' "2,2 + 4a2,3 

and have used the constant curvature condition aiRo — 2ao in the final line of (4.7). 

The final expression in (4.7) is a quadratic action with a Gauss-Bonnet term, a Weyl- 
squared term and a Ricci-scalar-squared term, together with a cosmological constant. The 
Gauss-Bonnet term can be dropped classically as a total divergence, leaving the action in 
the same quadratic form as discussed in our previous paper [1]. Reducing the action to a 
second-order form then follows exactly as in the case for f{R^^) actions, so that we obtain 
the same transformed action (3.9). Again we can verify that the transformed action has a 
vacuum solution at i? = i?o as is required for the consistency of our expansion. We also note 
that here too the expansion could have been made in terms of the variable Xx^j^^p, which 
could then be decomposed into its scalar and spin-two parts, the form of the decomposition 
depending on which vacuum is being considered. 

In conclusion, general actions of the form f{Rx^ivp) may have a variety of vacuum solutions, 
generically not apparent in a linear analysis. Around any vacuum the new degrees of freedom 
in the theory, aside from the massless graviton, can always be separated into a scalar field 
and a spin-two field. Unfortunately, the spin-two field is always ghost-like. 
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5. Higher-order Actions 



In this section, we shall briefly discuss how to extend our analysis to actions with higher- 
order equations of motion. We have seen that an action involving any function of the 
curvature tensors gives at most fourth-order equations of motion. For higher-order equations 
we need to consider actions which include some derivatives of the curvature. Here, we shall 
be concerned with only the simplest form. 



where we require that the function / has been reduced, by integration by parts and dropping 
total derivatives, so as to minimize k. The equations of motion following from similar actions 
have been considered by Schmidt [8] and Wands [9], who showed them to be equivalent to 
those for a set of scalar fields coupled to gravity. Here will shall show the full equivalence in 
a new way, working at the level of the action . 

It will be important to distinguish between possible forms of the function /. If we write 
/ = /(Ai, A2, . . . , Afc+i), (where in the action we have Ai = R, X2 = V^i? and so on), we find 
there are two possible cases: either df/d\k+i is a function of A^+i (case one) or it is not 
(case two). If it is not then it must be a function of A^, since otherwise the original form was 
reducible; that is, the action was not written in a form minimized with respect to k. Thus, 
in case two, we can always decompose / as 

case 2: /(Ai, A2, . . . , A^, A^+i) = 5'(Ai, A2, . . . , Xk)h+i + ^(Ai, A2, . . . , A^). (5.2) 

Further, we see that for case one the equations of motion are (4A; -|- 4)th-order, while for case 
two they are {4k + 2)th-ordcr. 

We would like to reduce the action (5.1) to a canonical second-order form, by introducing 
auxiliary fields. Given the order of the higher-derivative equations of motion we see that 
in case one we must introduce 2k + 1 new fields, while in case two we need only 2k new 
fields. The procedure we shall use is essentially a generahzation of Ostrogradski's method 
for reducing a higher-order action to a first-order form [12, 13], only that here we shall be 
reducing to a second-order form. (The Ostrogradski result is usual given as a Hamiltonian, 
but this can always be rewritten as a Helmhotz Lagrangian, the analog of the form we shall 
use.) 

We start by introducing a set of Lagrange multipliers, so 



Clearly eliminating fields via the Lagrange multiplier equations of motion, starting with /j.^ 
and working down in k returns one to the original action. 




(5.1) 




(V'Ai - A2) + • • • + A^fc (V'Afc - Afc+i)] . (5.3) 
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However, we have introduced at least one too many new fields. It is clear that the A^.+i 
equation of motion is purely algebraic and eliminating it will not introduce higher- derivatives; 
that is, the action will remain second order. The A^+i equation of motion reads 

{ f'^ (Ai, A2, . . . , Afc, Afc+i) in case 1, 
CAfc+i (5.4) 
5'(Ai, A2, . . . , Afc) in case 2, 

where we have distinguished between the two cases discussed above, and substituted the 
special form of / in case two. In case one, to form the analog of Ostrogradski's La- 
grangian, we solve the equation to give Afc+i as a function of Ai,...,Afc and fik, writing 
Afc+i = Afc+i(Ai, . . . , Xki^Jik)- It should be noted that in general the solution is not unique, 
and we must divide the original theory into pieces corresponding to different branches of the 
solution, just as in the case of actions of the form f{R) discussed in the previous section. In 
case two we simply substitute for /x^ and the special form of /. We get in case one 

case 1: S ^ ^ J d'^xy/^ /(Ai, A2, . . . , A^, Afc+i(Ai, . . . , A^, //fc)) + (i? - Ai) 

+^1 (V^Ai - A2) H h/xfc (^V^Afc- Afe+i(Ai,... ,Afc,Aifc)) , (5.5) 

the exact analog of the Ostrogradski-Helmhotz Lagrangian, while for case two we have a 
slightly different form, 

2: ^=^J c?'^V=^[MAi,A2,... ,Afc)+/.(i?-Ai) 

(V'Ai - A2) + • • • + ^(Ai, A2, . . . , Afc) V' Afc] . (5.6) 

We see that, as expected, in case one we have a total of 2A; + 1 auxiliary fields, while in case 
two we have only 2k new fields. 

All that remains is to transform the action into a canonical form, growing canonical kinetic 
energy terms for all the new auxiliary fields. Let us concentrate on actions of case one. Terms 
of the form //jV^Ai are easy to deal with. Simply introducing a pair of new fields, Xi i^i 

by 

\ = Xi + i>i: l^i = Xi- A: (5.7) 

we have, in the action, 

= j d'xV^ [- (Vxif + Wif] ■ (5.8) 

We see that the Xi field has a canonical kinetic term, but the ipi field has the wrong sign; 
it is a ghost. This is characteristic of higher- order theories. In reducing the theory to 



case 
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second-order the new fields always enter as a pair of a ghost-like field with a ordinary field. 
If we define a potential function, 

V{Xi, ■ ■ ■ ,Xk;i>i, - ■ ■ ,i>k) = ywAi + H1X2 H h /Ufc_iAfc + /ifeAfe+i(Ai, ... , Afc, Hk) 

~ /(-^i) -^2, • • • , Afe, Afc+i(Ai, ... , Afc, Hk)) (5.9) 

where it is understood that the right-hand side is evaluated at Aj = Xj + //j = Xj — ■0^, 
the action in case one becomes 



case 1: S — / d^x\p^ 



(5.10) 



To complete the transformation to canonical form all that is left is to make a conformal 
rescaling of the metric to remove the \xR coupling. As usual we define x = log \x and rescale 
to a metric g^^ = e^g^^, giving 

easel: 5=^ j d'x^ :r _ |(Vx)' - e'^^ ^ { (Vx.)' - (V^O'} 

-e-^^y(xi,...,Xfc;V'i,---,V'fc)]- (5.ii) 

We conclude that the original case one higher-derivative gravity theory is equivalent to 
canonical Einstein gravity coupled to 2A; + 1 scalar fields, A; + 1 of which, x and Xi, propagate 
physically and k of which, ipi^ are ghost-like. 

To put the case two action in canonical form is more complicated because of the 
g{\i, A2, . . . , Afc)V^Afc term. However, in principle, it is always possible to introduce a set 
of new fields {xi, • • • , Xk] i^i, ■ ■ ■ , V'fc} which simultaneously diagonalize the kinetic terms for 
the Aj and //j, though now the form of the transformation will depend on the function g. At 
least A; — 1 of the new fields will be ghost-hke. We can then make a conformal rescaling as in 
case one to put the action in the same canonical form (5.11), though the potential function 
will have a different form. 

The conclusion is that there is a procedure for rewriting the higher-order action (5.1) in 
a canonical second-order form. However, k of the new fields in case one and at least k — 1 
of the fields in case two will be ghost-likc. In each case we obtain a specific potential, and 
so we can again look for vacuum states as stationary points of the potential. Generically, as 
in the case of f{R) actions discussed in the previous section, the non-trivial vacua do not 
correspond to the flat-space solution of the original higher-derivative theory. 

As an simple example of this procedure consider the function f = R + a(V^i?)^. This 
example is of the flrst case since writing /(Ai, A2) = Ai + 0X2^ with \i = R and A2 = V^-R, 
we have df /d\2 = 2aX2, which is not independent of A2. Following our general procedure. 
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the first step is to introduce a set of Lagrange multipliers, 



1 

1 
2^ 



(fxsj—g 



J d^xyf^ [Ai + aX2^ + fx{R-Xi)+ m (V'Ai - As)] . (5.12) 

As discussed above we have introduced one too many auxiliary fields. We can eliminate A2 
by solving its equation of motion, which reads 

2a\2 - /ii = (5.13) 

implying A2 = ^i/2a. Substituting back into the action gives 

S = J d^x^/^ [^R + yUiV^Ai + Ai - jjiXi - \a^,i^] . (5.14) 

Next, to put the kinetic energy for Ai and jii in canonical form, we define Ai = xi + V'l and 

111 = Xi — i^i) so that 

3=^1 d'xV^g [liR - {Vxif + Wif - (/^ - 1) (Xi + V'l) - |« (Xi - ^if] ■ (5.15) 
Finally we make the conformal rescaling 'g^j^ = e^g^i, with % = log/i to put the action in 



canonical form 
1 



2fi:2 



d''x^ 



-9 



R - 1 (Vx)' - (Vxi)' + e-^ (Wi)' 



Thus we see that the original higher-derivative gravity theory is equivalent to canonical 
Einstein gravity coupled to three scalar fields. One field is ghost-like, and the potential has 
a single, unstable stationary point at % = Xi = = 0. 

As a second example consider / = a + fiR + ^fR^ + tKW'^R Writing Xi = R and A2 = V^-R, 
we now have df /dX2 = eAi, independent of A2, so this example clearly falls under case two. 
Repeating our procedure, we introduce Lagrange multipliers to give 

5 = -L y d^x^f^ [a + l3R + 7i?2 + eRV'^R] 

" 2^ / ^'^^^ + ^-^1 + ^-^i' + ^^^^^ + ii{R-Xi)+ 111 (V'Ai - A2)] . (5.17) 

Again we have too many new fields. The A2 equation of motion now reads 

eAi = 111, (5.18) 

which cannot be solved for A2 since the action is case two. However substituting this solution 
into the action eliminates both A2 and /xi, leaving the second-order form. 



S^^J d^x^ [liR - e (VAi)' - //Ai + a + /3Ai + 7A1'] . 



(5.19) 
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The kinetic energy for Ai is already in canonical form, but we must make a final conformal 
rescaling by = to put the action in the completely canonical form 

^=2^/ d'x^^-l(yxf -ee-^'i^Xif -V{\,,x)\ (5-20) 

where we have the potential, 

y(Ai, x) = e-^^ (e>^Ai - a - /3Ai - 7A1') (5.21) 

Again the original higher-derivative gravity is shown to be equivalent to ordinary Einstein 
gravity though now coupled to only two scalar fields. By choosing e > we can ensure that 
neither field is ghost-like. Note that we argued above that case two theories must have at 
least k — 1 ghost-hke fields. It is thus only in this special case of A; = 1 that we are able to 
have all the scalar fields non-ghost-hke. 

As before we obtain a specific potential for the fields. We find that ^(Ai, x) has a single 
stationary point at Ai = —2a/ (3, x — ^'^{P ~ '^<^l/P): provided P — Aa^/P > 0. Expanding 
around this point we find that, if a < and /3 < 0, we have a stable minimum. The 
value of the potential at the minimum is V^(Ai, x) = —<^/ {0^ — Aa^y) which is negative. We 
conclude that, for the given range of a, (3 and 7, the theory has a single stable vacuum 
state with negative cosmological constant. From Einstein's equation we find that this state 
is an anti-deSitter space with R = —2a/ — Aa^). Using the fact that R = e^R for 
covariantly constant we find that R = —2a/ (3. Thus the vacuum state also corresponds 
to an anti-deSitter space in the original higher-derivative theory. 



6. Conclusion 

The most important conclusion of this paper is that higher-derivative theories of gravita- 
tion generically have multiple stable vacua. One of these may be trivial, corresponding to 
fiat spacetime, but all the other vacua are non-trivial with the associated manifold being 
either deSitter or anti-deSitter spacetime with non- vanishing cosmological constant. While 
of interest from various points of view, such non-trivial vacua cannot represent the universe 
as it is now since the radius of curvature of these solutions is of order the inverse Planck 
mass. Thus, one might conclude that non-trivial gravitational vacua are irrelevant for par- 
ticle physics theories. However, this is not the case. We have recently shown that if we 
extend the methods and results of this paper to the realm of = 1 supergravity, then non- 
trivial vacua can exist with vanishing cosmological constant [14]. Furthermore, we find that 
supersymmetry is generically spontaneously broken in these vacuum states. It follows that 
higher-derivative A^ = 1 supergravitation could play a pivotal role in high energy physics. 
This possibility is being pursued elsewhere [15]. 
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